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Abstract 

Magnetic field influence on unsteady free convection flow of a second grade fluid near an infinite vertical flat plate with 
ramped wall temperature embedded in a porous medium is studied. It has been observed that magnitude of velocity as 
well as skin friction in case of ramped temperature is quite less than the isothermal temperature. Some special cases namely: 
(i) second grade fluid in the absence of magnetic field and porous medium and (ii) Newtonian fluid in the presence of 
magnetic field and porous medium, performing the same motion are obtained. Finally, the influence of various parameters 
is graphically shown. 
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Introduction 

The natural convection heat transfer from a vertical plate to a 
fluid has applications in many industrial processes. It was 
extensively studied by a number of researchers using different sets 
of thermal conditions at the bounding plate. Special mention can 
be made, for instance, to the studies of Raptis and Sing [1,2], 
Sacheti et al. [3], Chandran et al. [4,5] and Ganesan and Palani [6] 
that have determined analytical solutions for velocity and 
temperature using continuous and well-defmed conditions at the 
wall. Samiulhaq et al. [7] discussed the influence of radiation and 
porosity on the unsteady magnetohydrodynamic (MHD) flow past 
an infinite vertical oscillating plate with uniform heat flux in a 
porous medium. Keeping in mind the importance of shear stress 
on the boundary, Fetecau et al. [8] reinvestigated the problem of 
Samiulhaq et al. [7] by considering shear stress on the boundary. 
However, some practical problems may require non-uniform or 
arbitrary wall conditions. Chandran et al. [9] studied the unsteady 
free convection flow of an incompressible viscous fluid near a 
vertical plate with ramped wall temperature and compared the 
results with those of the plate with constant temperature. Recently, 
Seth and Ansari [10] and Seth et al. [1 1] found exact solutions for 
the MHD natural convection flow past an impulsively moving 
vertical plate with ramped wall temperature in the presence of 
thermal diffusion or radiation heat transfer. Narahari and Beg [12] 
considered the problem of Chandran et al. [9] for the impulsive 
motion of the plate with radiation and constant mass diffusion. 
More recently, Samiulhaq et al. [13] investigated the unsteady 
MHD flow past an impulsively started vertical plate present in a 
porous medium with thermal diffusion and ramped wall temper- 



ature. However, all aforementioned results refer to incompressible 
viscous fluids. 

Due to increasing significance of non-Newtonian fluids over the 
past few years, several researchers in the field are involved by 
valuable contributions in the study of flows of non-Newtonian 
fluids. It is due to their numerous applications in several areas, 
such as the plastic manufacture, performance of lubricants, food 
processing, or movement of biological fluids. These fluids are 
defined by a non-linear constitutive relationship between the stress 
and the rate of deformation tensors and, therefore, various models 
of non-Newtonian fluids have been proposed. Amongst them, the 
second grade fluids are the simplest subclass for which one can 
easily obtain analytical solutions. For these reasons and because, 
the second grade fluids can model many fluids such as dilute 
polymer solutions, slurry flows, industrial oils, many flow problems 
with various geometries and different mechanical and thermal 
boundary conditions have been studied. 

Tan and Masouka [14] investigated the Stokes' first problem for 
a second grade fluid in a porous half-space with a heated flat plate. 
Hayat and Abbas [15], by means of homotopy analysis method, 
have studied the heat transfer on the MHD flow of second grade 
fluids in a channel with porous medium. Closed form solutions for 
MHD flow of a second grade fluid through porous space are 
obtained by Khan et al. [16]. Thermal effects in Stokes' second 
problem for second grade fluid through a porous medium under 
the effect of magnetic field have been investigated by Srinivasa 
Rao et al [17]. Mustafa et al [18] have studied free convection flow 
of a viscoelastic second grade fluid along a vertical plate with 
power- law surface temperature. 

The influence of magnetic field is observed in several natural 
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and human-made flows. Magnetic fields are commonly applied in 
industry to pump, heat, levitate and stir liquid metals. There is the 
terrestrial magnetic field which is maintained by fluid flow in the 
earth's core, the solar magnetic field which originates sunspots and 
solar flares, and the galactic magnetic field which is thought to 
control the configuration of stars fi:'om interstellar clouds [19]. 
Three major technological innovations namely, (i) fast-breeder 
reactors used liquid sodium as a coolant which requires pumping; 
(ii) controlled thermonuclear fusion needs that the hot plasma be 
confined away from material surfaces by magnetic forces; and (iii) 
MHD power generation, in which ionized gas is propelled through 
a magnetic field were made by incorporating MHD in the field of 
engineering. The phenomenon concerning heat and mass transfer 
with MHD flow is important due to its numerous applications in 
science and technology Hayat et al. [20,21] and Hayat and Qasim 
[22]. The particular applications are found in buoyancy induced 
flows in the atmosphere, in bodies of water and quasi-solid bodies 
such as earth. Therefore, heat and mass transfer with MHD flow 
has been a subject of concern of several researchers (see for 
example, Katagiri [23], Jana et al. [24], Mandal and Mandal [25], 
Gosh [26], Jha and Apere [27], and the references therein). Other 
interesting results regarding the second grade fluids can be found 
in the references [28-34]. 

The purpose of this note is to extend some of the previous 
results to a larger class of fluids, namely to second grade fluids. 
More exactly, we establish exact solutions for velocity and 
temperature corresponding to the natural convection flow of a 
second grade fluid near an infinite vertical plate with ramped wall 
temperature. Apart from several other applications, the present 
study is significant and worthwhile as the exact solutions obtained 
in this paper are important not only that these solutions are new 
but as they can be used as checks for many approximate solutions 
and as tests for verifying numerical schemes. These solutions, 
obtained both for Pr 7^ 1 and Pr = 1, satisfy all imposed initial and 
boundary conditions. For comparison, the solutions corresponding 
to the plate with constant temperature are also established. Finally, 
temporal and spatial variations of velocity as well as those of the 
wall skin friction are graphically discussed. 

Mathematical Formulation of the Problem 

Let us consider the unsteady MHD flow of an incompressible 
second grade fluid near an infinite vertical plate with ramped wall 
temperature. The flow of electrically conducting fluid is taken in a 
porous medium. The x — axis is taken along the plate in the 
upward direction and y — axis is taken normal to the plane of the 
plate. A uniform magnetic field of strength Bq is acting in 
transverse direction to the flow as shown in Figure 1. Initially, at 
time / = 0, both the fluid and the plate are at rest to a constant 
temperature Too- At time / = 0^, the temperature of the plate is 
raised or lowered to Too -\-(T^; — Too)t/to when t<to, and 
thereafter, for t>to, is maintained at the constant temperature 
T^v The main purpose here is to study the free convection flow 
resulting from the ramped temperature profile of the bounding 
plate. 

It is assumed that the effects of viscous dissipation are negligible 
in the energy equation. One of the body force term corresponding 
to an MHD flow is the Lorentz force J x B. Where B is the total 
magnetic field and J is the current density. By using Ohm's law, 
the current density is given as 

J = a (E + V X B), (1) 

where cf is electrical conductivity of the fluid, E is the electric field, V 



is the velocity vector field, B = Bo+biwith Bq is the imposed 
magnetic field and biis the induced magnetic field. The current 
density J with the assumptions E = 0, bi=0 and 5 = 5o = (0,5o,0), 
where Bq is the strength of applied magnetic field Bq, modifies to 
J = (7 (V X Bo). Finally the Lorentz force becomes 
J X B = — cr^Q V,as mentioned by Hayat et al. [28] . For the problem 
under consideration, we assume the velocity of the following form 
\ = u(y,t)i where i is unit vector along x — axis. Under the usual 
Boussinesq's approximation of temperature gradient the equations 
governing the flow are: 

du(y,t) ( Oil d\ d^u(y,t) _ , 

(2) 

aBl ^ ^ <j> f ai d\ ^ ^ 

dT(y,t) ^ k d^T(y,t) 
dt pCp dy^ ' 

Here wis the velocity of the fluid in the x — direction, Tis its 
temperature, p is the density, gis the acceleration due to gravity, 
is the volumetric coefficient of thermal expansion, v is the 
kinematic viscosity, (j) is the porosity of the porous medium, ki is 
the permeability, k is the thermal conductivity, Cp is the specific 
heat of the fluid at constant pressure, and ai is one of the material 
module of second grade fluids. 

The initial and boundary conditions are: 

u(y,0) = 0,T(y,0) = To,; y>0, 

u(0, o = o,/>o, 

T(0,t) = Too^(T,,-Too)^ {ovO<t<to, (4) 

r(0,0 = ^w for/>/o, 

u(y,t)^0, T(y,t)^Tao Sisy^co Sindt>0. 



X 




Figure 1. Physical system and coordinate axes. 

doi:10.1371/joumal.pone.0088766.g001 
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0.04 




Figure 2. Comparison of velocity u(y,t) in Eq. (37) with Eq. (18) 
in Chandran et al. [9]. 

doi:1 0.1 371 /journal. pone.0088766.g002 

Introducing the following non-dimensional physical quantities 



y ^ t , 
y = r- ^ t = —,u 



^^,Pr = ^, (5) 




into Eqs. (2) and (3) and dropping out the notation we get 



Figure 3. Velocity profiles for different values of a with 

Pr =0.71, K = 2, M= 1 and ^ = 0.8. 
doi:10.1371/journal.pone.0088766.g003 

boundary conditions (8) will be determined by means of Laplace 
transforms. For comparison, the solutions corresponding to an 
isothermal plate with constant temperature are also established. 
Applying the Laplace transform to Eqs. (6), (7) and (8), we obtain 
the transformed equations 



d^u(y,t) d^u(y,t) du(y,t) , , , , ^ ^ 



dy^dt 



dt 



(1 + ag) 



^ ujy^q) 

dy^ 



-{cq + b)u{y,q) + e(y,q) = 0. 



(9) 



dt 



(7) 



-q^xQ(y,q) = ^; y,q>Q, 



(10) 



The adequate initial and boundary conditions are 



where t/(j,^)and 0(j,g)are Laplace transforms of w(j,^)and 
Q(y,i)^ together with the initial and boundary conditions in the 
transformed domain 



^(0,0 = 



tH{i)-{t-X)H(t-X) 



t/0;,0) = 0,%,0) = 0;j>0. 
w(0,0 = 0,/>0 

t , 0</<l 
1 , t>\~ 
wO;,0^0, %,0^0 asj^oo, />0, 

where H( )is the Heaviside step function and 



a= — , /? = M^ + — , c=l + — , M^= — ^to, 
fito K K p 



1 _ ve 



■to, to '- 



(8) 
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Solution of the Problem 

In the following, exact analytical solutions for the coupled 
partial differential equations (6) and (7) with the initial and 



-0.02': 



1 



Figure 4. Velocity profiles for different values of M with a = 0.4, 
Pr =0.71, K= \ and / = 0.8. 
doi:10.1371/joumal.pone.0088766.g004 
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whose inverse Laplace transform is given by 



u(0,q) = 0, 0(0,q)-- 



(11) 



i/(y,g)^0, 0(y,g)^O as y^co. 



The equation (10) is uncoupled to Eq. (9) and its solution with the 
corresponding conditions (11) is 

~e{ys) = {l-e-^)\\^^e-y^^^\ (12) 



d(y,t)=-e ^ 
a 



erfc (^^^ e-Ho ^(c-ab)zt^ dz 



erfc 



0 0 



f J e ^ lo(- \/(c — ab)zs] dsdz, 
\l^z) \ci J 



(19) 



and get 



Denoting by 



1 -JV^V^ 



' t erfc \- \ ^ e 



1/2(3^,0 = ^(^,0- erfc (i^^j;), 



(20) 



where Iq (g) is modified Bessel function of the first kind of order 
zero and erfc(g)is complementary error function. Consequently, 
(13) the expression for velocity in the (3;,0 — domain, can be written in 
the simple form 



and using the second shift property 
L-\e-""F{q)}=f{t-a)H(t-a) if f(t) = L-\F{q)), 



u(y,t)=U(y, t)H{t) -U{y,t-\)H(t-\). 



(21) 



where 



we obtain the following known result for the temperature 
distribution [9, Eq. (11)] 



U{y,i) = {u\ ® U2){i) = \u\{t — s) U2{y,s) ds, 

0 



e(y,t) = ei(y,t)-ei(y,t-\)H(t-iy 



(14) 



by 



The solution corresponding to Eqs. (9), (ll)i and (11)3 is given 



where the symbol {ui^U2){t) denotes convolution of Miand U2. A 
simple analysis clearly shows that both solutions (14) and (21), 
satisfy all imposed initial and boundary conditions. In order to 
show that w(0,0 = 0, for instance, we need the following integrals. 



u(y,q) = (l-e-^)UMU2(y,q), 



(15) 



where 



Uiiq)-- 



11 m2 



(xVxm2q{q + mxf-m2^ 



and 



U2{y,q)-- 



1 



(16) 



Pr-c , J(Fv-cf^4abFv 

with mi = ^ _ and m? = ^-^r 

2a Pr 2a Pr 

Laplace transform Ui{t) of Ui{q)is given by. 



The inverse 



^iW = ^[^1 smh(m2t) + m2COsh(m2t)]e-'^i^ - I. (17) 

In order to determine the inverse Laplace transform U2(y,t) of the 
function U2(y,q), we consider the following function: 



cq + b 
ag+ 1 



(18) 



0.15 



-0.05 



■ ■ ■ Isothermal temperature 
Ramped temperature 




1 



Figure 5. Velocity profiles for different values of K with a = 0.4, 
Pr =0.71,M = 0.5 and ^ = 0.8. 
doi:10.1371/journal.pone.0088766.g005 
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Figure 6. Velocity profiles for different values of Pr with a = 0.5, 

K = 2, M=l and / = 0.8. 

doi:10.1371/journal.pone.0088766.g006 
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-e ^ 
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e ^^\() ^- yj{c — aib)z dz = 

- a/(c — (xb)z ds dz = ^ (^l—e~c^^ . 



cz + s ^ 

e « lo 



Solutions for Plate with Constant Temperature 

In order to bring to light the effects of ramped temperature of 
the plate on the fluid flow, we must compare our results with those 
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Figure 7. Velocity profiles for different values of t with a = 0.5, 
Pr =0.71, M=l and K = 2. 
doi:10.1371/journal.pone.0088766.g007 



0.3 



o 



00 



1 


1 


I 


1 ^ 


M = 3 












-A. -A. .A. 

V V V 


a=0.4 - 








a=0.8 


^ 1 


1 


1 


a= 1.5 
1 


0 1 


2 


t 3 


4 t 



0.6 



^ 0.4 



CO 0.2 



1 1 

a = 0.75 


1 1 


o "1 n n o 

—A — A — A — A — A — * — ^ 




«^ M = 0.3 




eee M= 1.5 


1 1 


M = 2.5 
[ 1 



t 



10 



Figure 8. Variations of the skin friction for different values of a 
and M with K = l and Pr =1.5. 
doi:10.1371/journal.pone.0088766.g008 

corresponding to the flow near a plate with constant temperature. 
In this case, the initial and boundary conditions are the same 
excepting Eq. (8) that becomes ^(0,0= Ifor />0. The expression 
for the dimensionless temperature Q{y,i) is again the same 
obtained by Chandran et al. [9, Eq. (19)], i.e. 



%,0 = erfc( \^-\y 



(22) 



Introducing the expression of 0(y,O into Eq. (9), and following 
the same way as before we fmd that 



where 



1 



aPrg2 + (Pr -c)q-b' 
and its inverse Laplace transform is 



(23) 



(24) 
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(29) 



0.08 



Vdocily given by Eq. (2 1) 
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Figure 9. Comparison between values obtained with (21) and 
(38) for a=1.5, K = 2, M = 3, Pr =1.5. 
doi:1 0.1 371 /journal.pone.0088766.g009 



As regards the skin friction, in dimensionless form, is 

T^{t) = T(y,t)\y=o, 
where the shear stress T(j,/)is given by [28] 



(30) 



dtj dy 



Using equations (21) and (26) into above equation, we obtain: 
the shear stress for ramped temperature as 



x(y,t) = F(y,t)H{t) - F(y,t - \)H{t - 1), (31) 



where 



u^(t)=^ — sinh(m20^"'^i'. (25) ^ 

aPr m2 

Consequently, the dimensionless velocity corresponding to this The shear stress for constant temperature is 

case is 



u{yA = {mm2\t)= ^u,{t-s)u2{y.s)ds. (26) t(j,0-1^2(^-^)^^^^ (32) 



0 

where 



Nusselt Number and Skin Friction 

From velocity and temperature fields, the expressions for 
Nusselt number and skin friction can easily be determined. They 
are measures of the heat transfer rate and shear stress at the 
boundary. The Nusselt number Nu, can be written as 



Fx{t) = ux{t)+ sinh(m20^"'"i^ , 

Pr m2 

F2{t)= f^^^sinh(m20+ ^cosh(m20V""^'' 
\am2Pr Pr J 

and 
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Table 1. Absolute errors of velocity calculated by Eqs. (21) and (38). 



y 


Absolute errors 




|w(y,t)-v(y,t)| 




t = 0.5 


t=l.5 


t = 2.5 


t = 3.5 


0 


3.646 X 10"7 


1.315 X 10-6 


7.155 X 10-7 


5.379 X 10-6 


0.1 


2.648 X 10" 7 


6.287 X 10-7 


8.56 X 10-6 


3.056 X 10-5 


0.2 


4.775 X 10"7 


1.564x 10-6 


4.898 X 10-6 


2.116 X 10-5 


0.3 


3.346 X 10-7 


2.421 X 10-5 


1.791 X 10-4 


1.445 X 10-3 


0.4 


7.325 X 10"8 


6.559 X 10-6 


9.598 X 10-8 


1.461 X 10-5 


0.5 


3.41 X 10~7 


6.016X 10-6 


2.052 X 10-5 


1.909 X 10-4 


0.6 


6.605 X 10"7 


1.01 X 10-5 


4.413 X 10-6 


3.941 X 10-5 


0.7 


8.34x 10"7 


1.033 X 10-5 


7.164x 10-6 


5.199X 10-5 


0.8 


9.778 X 10"7 


1.181 X 10-5 


1.76 X 10-6 


4.563 X 10-6 


0.9 


1.016 X 10"6 


1.014x 10-5 


7.615 X 10-6 


5.402 X 10-5 


1.0 


1.005 X 10"6 


9.482 X 10-6 


4.201 X 10-6 


1.223 X 10-5 


1.1 


9.665 X 10"7 


8.622 X 10-6 


1.386 X 10-6 


2.128 X 10-5 


1.2 


8.803 X 10"7 


6.831 X 10-6 


3.171 X 10-6 


1.671 X 10-5 


1 ^ 


7 090 V 1 0- 7 
/ .yZy X iU / 


A. Q7^ V 1 n- 


Q 9 V 1 A- ^; 

y.z. X lU 0 


D.OUZ X lU J 


1.4 


7.055x 10-7 


4.327 X 10-6 


1.45x 10-6 


6.789x 10-6 


1.5 


6.17x10-7 


3.252X 10-6 


5.441 X 10-7 


2.456 X 10-6 


1.6 


5.32x10-7 


2.34x 10-6 


1.707X 10-8 


7.218X 10-7 


1.7 


4.5x 10-7 


1.515X 10-6 


2.137X 10-7 


3.462x 10-6 


1.8 


3.647x 10-7 


3.186X 10-7 


5.983 X 10-6 


4.18X 10-5 


1.9 


3.219x10-7 


1.368X 10-8 


5.342 X 10-6 


5.548x 10-5 


2.0 


2.817x10-7 


3.606 X 10-7 


6.712X 10-7 


1.315x10-5 



doi:l 0.1 371/journal.pone.0088766.t001 



Su2(y, t) 

Sy 



= \lvr w^''^' 1 '° 6 ^^^^^ )^^' 

0 



GO 

— ( 



lo ( - \/(c — cib)z s I ds dz. 



Ui(q)- 



1 



1 



m\=m2 = m-- 



aPr q^{q-\-m) 
Pr-1 
2a Pr 



/f 



and Pr 7^1, 



(33) 



Special Cases 

Te solutions corresponding to the flow of a second grade fluid 
with ramped wall temperature or constant temperature on the 
boundary in the absence of magnetic or porous effects can be 
immediately obtained from the general solutions (21) and (26) by 
making M^^O or K^co, respectively. However, if ^^00, the 
constant c = 1 and the corresponding solutions are different for 
Pr=land Pr^^l. So, for completion, we also give the exact 
solutions for velocity in two special cases. 

1. The case M = 0 and K^co 

By making M = Oand K^co, it results c=land b^O. The 
function U\{q) from Eq. (16)i becomes 



(34) 



The corresponding velocity u(y,t), after lengthy but straightfor- 
ward computations, is found to be 



ds 



(35) 



+- 



[m(r-^)-(l-£^-^^^-^>)]e-^erfc(^^^ h(^V^ dsdz, 



for PrT^ 1 respectively. 
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u(y,t)-- 



2a 



r 

60 6ocJ 
1 



(36) 



+ J(?-5)^e-^erfc(^^jlo(^^\/^^t/5^/z, for Pr = l. 

0 0 

2. The case a = 0 (Newtonian fluid with MHD and Porosity) 
For Pr = 1 



both cases of the plate heating. The fluid velocity is an increasing 
function of time t in the boundary layer then, finally it approaches 
to zero. In Fig. 8 are plotted the diagrams of the skin friction T^^{t) 
given by Eq. (31), for several values of the second grade parameter 
QL and magnetic parameter M. If both parameters increase, then 
the skin friction decreases. For a short time-interval the skin 
friction increases then approaches to a constant value. 

In order to obtain the closed form (21) of solution, we have used 
the Laplace transform method. In many problems, the inversion of 
image - functions can be a diflicult problem. Even if in our work, 
the inversion of function (15) is not too difficult, we present a 
numerical technique for inversion, namely, the Stehfest's algo- 
rithm [35]. Based on the Stehfest results, the inverse Laplace of the 
function H(y, q) is given by 



u{y,t) = ^ [a{y,t)H{t) -a{y,t-\)H{t-\)] (37) 



where 



hiy,t) = L-'{Hiy,q)}=^Y.djn{yJ^^) (38) 
where p is a positive integer and 



a{y,t)-- 



y+Oerfc 



In Figure 2, by making b^^it is observed that the graph of 
w(y,/)in Eq. (37) is similar to that of Eq. (18) from Chandran et al. 
[9] given by 



/bt 



\Ab 2 



)e^^erfch 



y 



-Vbt 



u{y,t) = uo{y,t)-uo{y,t-l)H{t-l), 



min(/» 



iP(2i)! 



i\(p-i)\(i-m-mi-jy 



(39) 



Here [r] denotes the integer part of the real number r and 



min(/»= If-PD 



(40) 



uo{y,t) 



where 



2t exp 



y t y^ 



At 



Numerical Results and Discussion 

The effects of different flow parameters have been analyzed by 
numerical calculations and graphical illustrations. A numerical 
algorithm was used in order to compare the analytical solutions 
with the numerical solutions. 

The velocity field, for various values of second grade parameter 
cx, is described in Fig. 3. The effect of the second grade parameter 
is to decrease velocity throughout the flow field when cx increases. 
It is also clear that, the velocity approaches to zero at the far away 
from the plate. It is noticed that, the thickness of the boundary 
layer increases if the second grade parameter decreases. For 
ramped temperature on the plate, fluids flow slower than for the 
constant plate temperature. The effect of the magnetic strength on 
the motion of the fluid, for both heating cases, is analyzed in Fig. 4. 
Increasing of the magnetic parameter decelerates the motion of 
the fluid in the boundary layer. Therefore, the magnetic field acts 
like a drag force. The influence of the permeability parameter K is 
shown in Fig. 5. It is observed that the velocity field is an 
increasing function of K. As expected, the increase of the 
permeability of the porous medium reduces the drag force and, 
therefore, fluid velocity increases. The effect of Prandtl number on 
the velocity field is sketched in Fig. 6. It is also clear that, the 
increase of the Prandtl number decelerates the motion of the fluid. 
In Fig. 7 are plotted the diagrams of velocity u(y, t), versus t, for 



Applying the formula (38) to Laplace transform u{y,q) given by 
equation (15), the values of velocity field u(y, t) are obtained. As 
shown in Figure 9 and in Table 1, the values ??of the function u(y, 
t) obtained by formulae (21) and (38) are in excellent agreement. In 
Table 1 we denoted by u(y, t) and v(y, t) the values of velocity given 
by Eq. (21), respectively, by Eq. (38). The Table 1 contains the 
absolute errors \u{y ,t) — v{y ,t)\. Figure 9 confirms the velocity 
properties shown in Figure 7, namely the increasing of velocity in 
the boundary layer when, the time increases. 

Limitations of the Study and Future 
Recommendations 

It is important to bring to light various limitations of this 
research. A discussion of these limitations will not only assist 
readers to understand this study, but also provide an opportunity 
to extend the current research. The following assumptions and 
limitations are considered 

• Flow is incompressible and laminar. 

• Flow is one dimensional and uni-directional. 

• A uniform magnetic field is applied outward direction 
perpendicular to the flow. 

• It is assumed that the effects of viscous dissipation in the energy 
equation are negligible. 

• Electric field due to polarization of charges is not considered. 

The mathematical model of second grade fluids offers, in 
general, possibilities to find of analytic solutions. Unfortunately, 
this model does not exhibit some significant features of some fluids. 
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From this reason the present work can be extended to other more 
complex models, such as the power - law fluids of second grade in 
which the fluid may exhibit normal stresses, shear thinning or 
shear thickening behavior. Also, the approaches of some models 
with fractional derivatives in various geometrical configurations 
can be interesting. The present study provides analytical solutions 
in the closed form which can be used as a benchmark by 
numerical analysts. 

Conclusions 

Exact solutions corresponding to the ramped wall temperature 
of unsteady MHD free convection flow of a second grade fluid in a 
porous medium are established. Solutions are obtained by using 
Laplace transform technique. The obtained solutions can easily be 
reduced to similar solutions for Newtonian fluids. They can be 
used to develop new exact solutions corresponding to free 
convection flows of several non-Newtonian fluids. The corre- 
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